Radiation pressure on a moving body: beyond the Doppler effect 
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The dependence of macroscopic radiation pressure on the velocity of the object being pushed is 
commonly attributed to the Doppler effect. This need not be the case, and here we highlight velocity 
dependent radiation pressure terms that have their origins in the mixing of s and p polarizations 
brought about by the Lorentz transformation between the lab and the material rest frame, rather 
than in the corresponding transformation of frequency and wavevector. The theory we develop is 
relevant to the nano-optomechanics of moving bodies. 
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The force of radiation pressure, which arises from the 
momentum of light, depends upon the velocity of the 
body being pushed. As an example, take light incident 
onto the surface of a perfectly reflecting mirror when this 
mirror is in motion in the direction of the surface normal. 
In the rest frame, the momentum and rate of arrival of 
the incident photons will appear either increased or de- 
creased relative to the lab, depending whether the mirror 
moves towards, or away from the light source. This is a 
consequence of the Doppler effect, and implies that the 
velocity dependent part of the radiation pressure acts to 
slow the mirror down. When attached to an oscillatory 
degree of freedom this effect appears as a friction-like 
term in the equation of motion and is therefore termed 
radiation damping [1] [5] . If the ideal mirror is replaced 
by a material that is highly dispersive (e.g. a Bragg mir- 
ror), the reflection coefficient of the moving medium is 
also strongly dependent upon the Doppler shifted fre- 
quency, and this can turn damping into heating [31 13] . 

Another instance of a velocity dependent radiation 
pressure effect is universal drag whereby all media in 
relative motion at finite temperature experience a fric- 
tional force due to the electromagnetic field [5]; this too 
has its origins in the Doppler effect. One consequence of 
universal drag is that when two polarizable bodies with 
equal rest frame temperature (T ^ 0) move relative to 
one another, they act to slow one another down via the 
electromagnetic field. 

The relative motion of a dielectric medium, however, 
has a peculiar effect on the polarization of light as 
well as its propagation. The moving medium behaves 
as a magnetoelectric |6l [7], which implies a change in 
the polarization composition (s & p) upon interacting 
with the medium. Here we raise the question whether 
there exist optical forces, at least from a classical stand- 
point 16 , arising from polarization mixing effects in a 
moving medium, in spite of the fact that many features 
of reflection from, and transmission through moving di- 




FIG. 1: In the lab frame the S matrix connects input 
fj') and output (ftqL iCtq^') amplitudes of fields scat- 
tering from the thin slab moving laterally with velocity V . 
The s and p polarization unit vectors, respectively, ei and e.2 
are here displayed for the special case of a plane wave prop- 
agating in the (x — y) plane with fc|||| y. More generally the 
unit vectors ei and e^-' lie in the (z — y) plane and in the 
incidence plane, i.e. the (a; — fe||) plane, respectively. 



electrics have been considered before ^i9^. What's more, 
is there an interplay between polarization mixing effects 
and Doppler effect that is important to the radiation 
pressure experienced by a moving dielectric? Here we 
show that this is indeed the case, and that there are sit- 
uations where polarization mixing effects occur at the 
same order of V/c as the Doppler effect, and hence cannot 
be neglected. Our theoretical framework, based on gen- 
eralized scattering matrix methods, is readily amenable 
to single out these non-Doppler contributions to the ra- 
diation pressure experienced by a moving dielectric. We 
further examine a clear-cut configuration where the effect 
of radiation pressure is solely due to polarization mixing. 



The lab frame electric field of radiation exerting pres- 
* Electronic address: |sarh@st-andrews . ac . uk| sure On a moving dielectric slab (see figure [ij can be 



2 



written as, 

E{x,t)^ g(±)^^±)gZ(±|fe.|a.+fc||.x-<^t). (1) 

±,9=1,2 

where the surface has a normal parallel to 5;, fc|| = 
kyjj + k^z and fc**' = ±|fca;|i + fey, with all wavevec- 
tor components real and with polarization unit vectors 

e^*' = Ci = i; X fc||, 62*' = k X ei, that are equiv- 
alent to s (g = 1) and p {q = 2) polarization, respec- 
tively. The electric field transformation into the rest 
(primed) frame of the moving slab depends upon its di- 
rection of motion. In particular, for {V = Vxx), the 
rest frame field amplitudes, a')*' are found to be re- 
lated to the lab frame by a constant of proportional- 
ity, a';±' = (c^"±Vc.)af', with t^'<±' = 7(w T 
and 7 = (1 — V^/c^)~^/^. Meanwhile for lateral motion 
{V = Vyi)) we find, a' — {uj' /uj)M • a, where, 

,' ( }l ^'A, (2) 

Uj' = 7(0; - Vyky), rj = \kx\kz/{k1^ - VyUlky/c'^), Oi^ = 

(a^j^^', aj^"', aJj^', Qfj'), and cTz is the usual Pauli matrix. 
Note that all of the particular velocity dependent radia- 
tion pressure effects mentioned in the introduction arise 
from the Doppler shift, and this holds also here for mo- 
tion along X, where the field amplitudes are only mod- 
ified by the factors w''*'/'^- The situation changes for 
dielectric media in lateral motion where the amplitudes 
are modified both by the Doppler shifted frequencies and 
by the transformation The unitary 4x4 matrix M 
describes how the two polarizations are coupled or mixed 
by the motion of the medium, and this specifically hinges 
on the off diagonal 2x2 matrices within ([2| that depend 
on the lateral velocity Vy and the parameter rj. Such a 
polarization mixing effect vanishes for fc^ = in which 
case the incidence plane contains the y axis along which 
the slab moves and thus it remains a plane of mirror 
symmetry for the system assuring that s and p polar- 
izations be independent. Conversely, when fc^ ^ and 
fcy = VyLoky/c^, we have rj — )■ cx), and the polarization is 
now completely converted from one type to another. We 
wish to investigate to what extent mixing of the s and p 
polarization composition affects radiation pressure when 
light is incident onto a moving medium. We will then 
consider in what follows only the case of a laterally mov- 
ing slab of thickness, d, centred at x,y,z = 0, and of y-z 
area, A = LyL^, where Ly,Lz ^ d. In both free space 
regions, x < (L), and x > d/2 (R), we assume an 

electromagnetic field in the form of a right-going wave 
plus a left-going wave as in ([l]). In the case of a slab in 
lateral motion, the frequency is conserved upon refiection 
and transmission, and so we consider all waves to have a 
single frequency, uj. The average four-force experienced 
by the slab can be written as the difference between in- 
put and output photon flux in each polarization, directed 
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FIG. 2: (a) Normal force (|9| as a function of the angle of 
incidence 6, acting on a slab of thickness 5c/cj with (a) e = 
5.0 + 0.011 and (b) e = 0.001-l-O.OOli and fj. = 1. The incident 
wave-vector is written as, k = fea;a; -I- fc||fc|| — (to / c)[cos{9)x + 
sin(6')(sin(x)y -I- cos(x)z)]. Inset shows the s (dashed) and p 
(solid) reflection coeflicients for such a slab . 



along the wave-vector, k'^ — (w/c, ij/cj,!, fey). 




This expression can be derived from a consideration of 
the energy-momentum tensor averaged over a time inter- 
val, T 3> uj~^, as in [3.,. In general, we have knowledge 
of only the input flelds, a.^^ = (a^"^', a[~^, a2t!i '^2~r)i 
not the output flelds, Oi^„.^ = {a\^^, a^"' , aj"^, ckj'l )• The 
relationship between the two is given by the scattering 
matrix, S: ckqut = S • ckin- In our case the motion of 
the medium will mix the polarizations, and we must deal 
with a 4 X 4 scattering matrix, rather than the standard 
2x2 object ^U\, the elements of which are written in 



3 



general as, 



5^21 5^22 



(Til Till Ti2 ni2\ 

Ti-ii Til TZi2 Ti2 

T21 Ti-21 T22 T^22 

\Ti-2l T21 1^22 T22 ) 



(4) 



In the case of media which do not mix s and p polariza- 
tions, Q reduces to a direct sum, S = Sn ® 822- 

Using Q, the radiation pressure experienced by the 
moving slab can be written in terms of the incident field 
amplitudes and the rest frame scattering matrix, S'{k'J^). 
First, (H is written in the rest frame of the medium in 
terms of S' and the field amplitudes, o:',^. After perform- 
ing a Lorentz transformation back to the lab frame, and 
applying ([2]), we find the four-force in the lab frame in 
terms of lab frame field amplitudes. 



dt 



a! 



( [u - S^S) LJ/C \ 



{R-S^RS)\k^ 
(I4 - S^S)ky 
V {l4-~S^S)k, 



a. 



(5) 



where, i? = cr^ © cTz , and the field amplitudes outside of 
the vector multiply the matrices inside. The lab frame 
scattering matrix, S is related to the rest frame ma- 
trix by, S = M^'S'M. The conclusion we can draw 
from comparing ([s]) with its limit when Vj^ — > is that 
a laterally moving medium responds akin to a station- 
ary medium with a peculiar frequency response, and re- 
flection and transmission coefficients that are related to 
the rest frame reflection and transmission coefficients by 
a unitary transformation. At a frequency w' where the 
material has negligible loss, S'^ S' — I4, all components 
of the four-force, besides (dP^/dt), vanish. Polariza- 
tions are not mixed in the rest frame where one has 

s' = s'li © s'22 with 7^^ = = r[, ku^n^ t[, 

where r[ ^ ^-iid t'-^ ^ respectively the rest frame re- 
fiection and transmission amplitudes for s and p polar- 
ization. In this case one has the lab frame scattering 
matrix. 



{S'li^T z — <y zS'2 




22) 



(6) 



where the off-diagonal blocks linear in {riVy/c) are re- 
sponsible for the polarization mixing effects. 

The four-force in ^ along with ^ is an important 
result, as it embeds a new radiation pressure effect stem- 
ming from mixing of the s and p polarizations associated 
with moving frames. Two specific examples are now il- 



lustrated. First, take the case of a plane polarized wave 
propagating in the direction k onto the left of a later- 
ally moving slab, the two input polarization amplitudes 
being a^^' and a'"2l- The resulting normalized [17] nor- 
mal [x) force is. 



'd_n 

, dt 



l + ai 

W2 



ai 



VVy 



-a2 



0.2 



0L2 



VVy 



-ai 



{2R[+A[)\ai\^ + {2R'2+A'2)\a2\ 



\kx\kz Vy 



{2Ri + Ai - 2R2 - A2){ala2 + c.c.) (7) 



where we define ai — a^^' / y'lai^L P + |q^2^l — 

i\r'f-\t'f)/{l+vX/c^) = {2R[+A[~l)l{lWV^/^). 
The last expression is correct to first order in Vy/c 
(r] ~ \kx\kz/k?^). Other terms linear in Vy/c would ap- 



pear upon expressing R[ and A'^ in terms of the lab frame 
quantities through the transformation of frequency and 
wavector components. Such an expansion would bring 
out familiar Doppler terms that are independent of the 
polarization mixing terms proportional to {ala2 + c.c.) 
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shown on the right hand side of ([t]). lossless medium {A'^ = 0), takes a similar form, 

The normalized lateral (y) force, which vanishes for a 



. dt / 



c^\kx\ ky 



ai 



Wy 



0,2 



OL2 



c I kx I ky 



Wy 



A[\ai\^ + A'2\a2\ 



{Ai - A2){ala2 + c.c.) (8) 



with 6, ^ (- \r'f - \t'f)/{l + rjX'/c') = (^'. " 1)/(1 + 
rf'Vy /c^). In general, polarization mixing contributions 



to the velocity dependent part of the force (7][8) are com- 
parable to the Doppler terms [18| . Since they are propor- 
tional to (0^0:2 + c.c), they vanish for unpolarized light, 
such as e.g. thermal radiation, where ai and a2 do not 
have a definite mutual phase relation, and also for cir- 
cularly polarized light {ai = l/-\/2, 012 = ±i/-\/2). Even 
when ?7 7^ 0, the polarization mixing effects vanish to lin- 
ear order in Vy/c if in the lab frame the light is cither 
purely s or purely p polarized, yet they survive and could 
be easily modulated or even flipped in sign by rotating 
the linear polarization when reflection and transmission 
for s and p polarizations are significantly different. In 
particular, for linearly polarized light with ai = sin (p 
and a2 = cos ip, corresponding to an electric field whose 
polarization vector makes an angle ip with e^^\ one has 
(a*a2 + c.c.) = sin(2(p). 

As a second example, we consider a situation where 
radiation pressure arises only when both the mixing 
coefficient t] and the velocity do not vanish. Take 
two linearly polarized plane waves of equal intensity 
propagating m the direction k and k , respec- 
tively onto the left and the right of a laterally moving 
slab. In much the same way as done above, the in- 
put field amplitudes may be chosen so that a^^ — 
f (e*^^ sin(i^L), e"'^ sin^ipn), e"-^ cos{(Pl), e"'"' cos{(Pr)), 
where AC = Cl ~ Cr is the phase mismatch between 
the two input waves. Note that with our choice of 
polarizations the plane wave coming from the right 
and having ipn = n — ip^ is the mirror image of that 
coming from the left with respect to the {z-y) plane of 
symmetry. Taking an average over as appropriate 
for left and right incoming waves not having a definite 
mutual phase relation, the normal component of the 
force in this case can be obtained from ([t]) and the use 
of the mirror symmetry just described. In particular, 
for the specific choice ipR — lpl = 'P and retaining only 
terms linear in Vy/c, one has 



dt 



2c^kx \kx\k2, Vy 



2(i?i-i?2)+Ai 



^2] sin(2v3) 
(9) 



This force, that is present even when equal incident 
intensities are used (in both polarization channels) on 
both sides of the slab, is entirely due to the asymme- 
try in the polarization transformation between left and 
right propagating waves as described by ([2]) and does not 
seem to have been recognised before. In addition, the 
force ([9|, which does not contain Doppler contributions, 
is typically of the same order as the velocity dependent 
radiation pressure stemming from the Doppler effect as 
mentioned in the introduction. Figure |2^ shows the de- 
pendence of this force on the angle of incidence for a thin 
dielectric slab when p = ■n/A. Again, there are clearly 
situations that make (|9| vanish altogether, such as when 
we have equal reflection and absorption for the two po- 
larizations, e.g. for surfaces behaving as ideal mirrors. 



For small (but not zero) fcy, as ry ~ |^a:|^z/fcjj becomes 
very large, there is a large degree of polarization mixing. 
For ordinary materials this does not have an effect due to 
the near degeneracy of s and p reflection and transmission 
coefScients at small angles of incidence. However, there 
are materials where the reflection and transmission coef- 
ficients remain significantly different close to normal in- 
cidence. For example, such a dependence on angle could 
apply to epsilon near zero metamaterials [111112] . and fig- 
ure shows more than an order of magnitude increase 
in the polarization mixing contribution to the velocity 
dependent force. 

The configuration leading to (|9| is of special interest as 
it enables one to observe pure polarization mixing contri- 
butions to the velocity dependence of radiation pressure. 
This, in turn, can simply be tuned via the polarization 
angle, for instance it can be turned on and off by setting 
PL = 7r/4 and switching p^ respectively to 7r/4 and 37r/4. 
Such an ability can be important in the emerging field of 
nano-optomechanics [131, where small optical forces can 
be used to manipulate the centre of mass of mechanical 
oscillators of very light masses, or even atomic clouds [3]- 
In particular, it can allow to control the dynamics of a 
two-dimensional oscillator in which the motion along one 
axis (y in our case) can be optomechanically coupled to 
that along the other one {x in our case). 
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